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Abstract 



In a previous work, we prove the existence of weak solutions to an initial-boundary 
value problem, with H 1 ^) initial data, for a system of partial differential equa- 
i-£h , tions, which consists of the equations of linear elasticity and a nonlinear, degenerate 

parabolic equation of second order. Assuming in this article the initial data is in 
H 2 (£l), we investigate the regularity of weak solutions that is difficult due to the 
gradient term which plays a role of a weight. The problem models the behavior in 
time of materials with martensitic phase transitions. This model with diffusive phase 
interfaces was derived from a model with sharp interfaces, whose evolution is driven 
by configurational forces, and can be thought to be a regularization of that model. 
■ Our proof, in which the difficulties are caused by the weight in the principle term, is 

0^ , only valid in one space dimension. 
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1 Introduction 



- 1—1 

X 

Many inhomogeneous systems can be characterized by domains of different phases sep- 
arated by a distinct interface [16]. Driven out of equilibrium, their dynamics result in 
the evolution of those interfaces which might develop into structures (compositional and 
structural inhomogeneities) with characteristic length scales at the nano-, micro- or meso- 
scale. To a large extent, the material properties of such systems are determined by those 
structures of small-scale. It is thus important to understand precisely the mechanisms 
that drive the evolution of those structures. In this article we are interested in a model 
for the evolution, driven by configurational forces, of microstructures in elastically de- 
formable solids. Materials microstructures may consist of spatially distributed phases of 
different compositions and/or crystal structures, grains of different orientations, domains 
of different structural variants, domains of different electrical or magnetic polarizations, 
and structural defects, see e.g. [15]. These structural features usually have an intermedi- 
ate mesoscopic length scale in the range of nanometers to microns. The size, shape, and 
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spatial arrangement of the local structural features in a microstructure play a critical 
role in determining the physical properties of a material. Because of the complex and 
nonlinear nature of microstructure evolution, numerical approaches are often employed. 

There are two main types of modeling for the evolution of microstructures. In the 
conventional approach, the regions separating the domains are treated as mathematically 
sharp interfaces. The local inter facial velocity is then determined as part of the boundary 
conditions, or is calculated from the driving force for interface motion and the interfacial 
mobility. This approach requires the explicit tracking of the interface positions. Such an 
interface-tracking approach can be successful in one-dimensional systems, however it will 
be impractical for complicated three-dimensional microstructures. Therefore, during the 
past decades, another approach has been invented, namely, the phase-field approach in 
which the interface is not of zero thickness, instead an interfacial region with thickness 
of certain order of a small regularization parameter. Though it is still a young discipline 
in condensed matter physics, this approach has emerged to be one of the most powerful 
methods for modeling the evolution of microstructures. It can be traced back the theory 
of diffuse-interface description, which is developed, independently, more than a century 
ago by van der Waals [27] and some half century ago by Cahn and Hilliard |13| . 

The two well-known models for temporal evolution of microstructures are the Cahn- 
Hilliard/ Allen- Cahn equations corresponding, respectively, to the case that the order 
parameter is conserved and not conserved. These phase field models describe microstruc- 
ture phenomena at the mesoscale, and one suitable limit of it may be the corresponding 
sharp- or thin-interface descriptions. In this article we study a model for the behavior in 
time of materials with diffusionless phase transitions. The model has diffusive interfaces 
and consists of the partial differential equations of linear elasticity coupled to a quasilin- 
ear, non-uniformly parabolic equation of second order that differs from the Allen-Calm 
equation (the Cahn-Hilliard equation in the case that the order parameter is conserved) 
by a gradient term. It is derived in [H [4] from a sharp interface model for diffusionless 
phase transitions and can be considered to be a regularization of that model. To ver- 
ify the validity of the new model, mathematical analysis has been carried out for the 
existence of weak solutions to initial boundary value problems in one space dimension, 
[U El [28] , the existence of spherically symmetric solutions [8] , the motion of interfaces 
[6], and the existence of traveling waves [21J. In the present article, the existence and 
regularity of weak solutions to an initial-boundary value problem will be studied. We 
first formulate this initial-boundary value problem in the three-dimensional case, reduce 
it to the one-dimensional case and prove the existence of weak solutions to this one di- 
mensional problem, then study the regularity of weak solutions by assuming that the 
initial data is in H 2 (£l). 

Let Q C M 3 be an open set. It represents the material points of a solid body. The 
different phases are characterized by the order parameter S(t, x) 6 R. A value of S(t, x) 
near to S- (which is a real number) indicates that the material is in the matrix phase at 
the point x € O at time t, a value near to S+ indicates that the material is in the second 
phase. The other unknowns are the displacement u(t, x) £ M 3 of the material point x at 
time t and the Cauchy stress tensor T(t, x) € S 3 , where <S 3 denotes the set of symmetric 
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3 x 3-matrices. The unknowns must satisfy the quasi-static equations 



div x T(t, x) 
T(t,x) 

S t (t,x) 




for (t, x) £ (0, oo) x $7. The boundary and initial conditions are 



u(t, x) 
S(0,x) 



7(t, x), S(t, x) = 0, (t, x) e [0, oo) x dtt, 
Sq(x), x £ O. 



(1.4) 
(1.5) 



Here V x u denotes the 3x3-matrix of first order derivatives of u, the deformation gradient, 
(S7 x u) T denotes the transposed matrix and 



is the strain tensor, e £ <S 3 is a given matrix, the misfit strain, and D : <S 3 — >■ 5 3 is a 
linear, symmetric, positive definite mapping, the elasticity tensor. In the free energy 



we choose for ip £ C 2 (M, [0, oo)) a double well potential with minima at S = S- and 
S = S + . The scalar product of two matrices is A ■ B = ^ aijbij. Also, ips is the partial 
derivative, c > is a constant and v is a small positive constant. Given are the volume 
force b : [0, oo) x O -> M 3 and the data 7 : [0, 00) x dQ -> M 3 , : O -> E. 

This completes the formulation of the initial-boundary value problem. The equations 
(II. ID and (jl.2p differ from the system of linear elasticity only by the term eS. The 
evolution equation fjl .3|) for the order parameter S is non-uniformly parabolic because of 
the term z/AS , |V a; S'|. Since this initial-boundary value problem is derived from a sharp 
interface model, to verify that it is indeed a diffusive interface model regularizing the 
sharp interface model, it must be shown that equations (jl-lj) - (11.51) with positive v have 
solutions which exist global in time and is more regular if the initial data is more regular, 
and that these solutions tend to solutions of the sharp interface model for v — > 0. This 
would also be a method to prove existence of solutions to the original sharp interface 
model. 

In this article we show that in one space dimension the initial-boundary value problem 
has solutions, and study the regularity of these weak solutions with H 2 (Q) initial data. 
Whether solutions in three space dimensions exist and whether these solutions converge 
to a solution of the sharp interface model for v — > is still an open problem to be 
investigated later. The model and therefore the existence result is of interest not only 
in three dimensions but also in one space dimension. Moreover we believe that this one- 
dimensional existence result can also be helpful for an existence proof for higher space 
dimensions. 

Related to our investigations is the model for diffusion dominated phase transforma- 
tions obtained by coupling the elasticity equations (jl.ip . (jl.2p with the Cahn-Hilliard 
equation. This model has recently been studied in 114 1 IB]- 





(1.6) 
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Statement of the main result. We now assume that all functions only depend on the 
variables x\ and t, and, to simplify the notation, denote x% by x. The set Q = (a, d) is a 
bounded open interval with constants a < d. We write Qt b '■= (0,T e ) x f2, where T e is a 
positive constant, and define 

(v,<p)z = / v(y)(p(y) dy , 
Jz 

for Z = Q or Z = Qt £ - If v is a function defined on Qj< e we denote the mapping 
x i — y v(t,x) by v(t). If no confusion is possible we sometimes drop the argument t and 
write v = v(t). We still allow that the material points can be displaced in three directions, 
hence u(t,x) G M 3 , T(t,x) G 5 3 and Sell. If we denote the first column of the matrix 
T(t,x) by Ti(t,x) and set 

e{u x ) = X - (K.,0, 0) + (u x , 0, 0) T ) G S 3 , 

then with these definitions the equations (jl.ip - (jl.3p in the case of one space dimension 
can be written in the form 

-T lx = b, (1.7) 
T = D(e(u x )-eS), (1.8) 

# = cfr-e-V>'(5) + ^) 15,1, (1.9) 



which must be satisfied in Qt £ - Here we have inserted if>s{e, S) = —T ■ e + tp'(S). Since 
the equations (jl.7p . (jl.8p are linear, the inhomogeneous Dirichlet boundary condition for 
u can be reduced in the standard way to the homogeneous condition. For simplicity we 
thus assume that 7 = 0. The initial and boundary conditions therefore are 

u(t, x) = 0, (t, x) G (0, T e ) x dn, (1.10) 
S(t,x) = 0, (t, x) G (0, T e ) x dn, (1.11) 
5(0, x) = So(x), xen. (1.12) 

To define weak solutions of this initial-boundary value problem we note that because 
of \{\y\y)' = \y\ equation (II. 9|) is equivalent to 

St-CLs\{\S x \S x ) x -c(T-e-j>'{S)) 15*1=0. (1.13) 

Definition 1.1 Let b G L°°(0, T e , L 2 (Q.)), 5 G L°°(S1). A function (u,T,S) with 

u G L°°(p,T e ;WQ ,00 (Q)) : (1.14) 
T G L°°(QT e ), (1.15) 
5GL oo (Q Te )nL oo (0,T e ; J ff 1 (f])), (1.16) 

is a u>ea£; solution to problem ji.7] ) - hi. 12(1 , if equations jl.?] ) - luzf/i hi. 10(1 are 

satisfied in the weak sense, and if for all <p G C^°((— 00, T e ) x O) 

(5, ^)Q Te "^^(15,15,, ^) QTe +c ((T ■ e - ft(S)) \S X \, <p) +(5 , <^(0)) n = 0. (1.17) 
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We make the following assumption on the nonlinearity ip(S). 
Assumption for ip(S): 

ip(S) G C°°(R) is a double- well potential which has two local minima at 
S- and S + with S_ < S + and one local maximum at 5* with S_ < S* < S + ,(1.18) 
and satisfies ip'(S) > for S_ < S < S* and if/(S)) < for S* < S < S+. 

One typical example is: ifj(S) = (5(1 — 5)) 2 with 5_ = 0, 5+ = 1. 
Now we are in a position to state the main result of this article. 

Theorem 1.2 Assume that $(S) satisfies (PTTT81) . Then for all 5 G H^(Sl) f] H 2 (Sl) and 
b G C(Q Te ) with bt G C(Q T J, there exists a weak solution (u, T, S) to problem fli. 7| j - 
111.12]) . which in addition to {1. - ( [i.J7| j satisfies 

UtEC^T^H 1 ^)), T t eC([0,T e ];L 2 (Q)) (1.19) 

and 

S t GL°°(0,r e ;L 2 (O)), 15,15, er^T^ 1 ^)), S x € L°°(Q Tlt ). (1.20) 

The remaining sections are devoted to the proof of this theorem. The main difficulty 
in the proof stems from the fact that the coefficient z^|5 x | of the highest order derivative 
S xx in equation (I1.9P is not bounded away from zero and that it is not differentiable with 
respect to S x . 

To prove Theorem 11.21 we therefore consider in Section 2 a modified initial-boundary 
value problem which consists of (|1.7p . (|1 .8[) . (jl.lOp - (|1 . 12[) and the equation 

S t -cv\S x \ K S xx -c(T-e-i>'(S?)(\S x \ K -K) = 0, x G SI, t > (1.21) 

with a constant k > 0. Here we use the notation 

[p|« : = V^+W- (1-22) 

Since (ll.2ip is a uniformly (for |5 X | < M, M is a positive constant) parabolic equation we 
can use a standard theorem to conclude that the modified initial-boundary value problem 
has a sufficiently smooth solution (u K ,T K , 5 K ). For this solution we derive in Section 3 a 
priori estimates independent of k G (0, 1]. 

The function \p\ is smoothed by \p\ K in (jl.22p which is different from that in [3|, thus 
we also prove the existence of weak solutions, though our main concern of this article is 
the regularity of solutions. To select a subsequence converging to a solution for n — > we 
need a compactness result. However, our a priori estimates of S xx depend on a weight 
\S X \ K , and are not strong enough to show that the sequence S£ is compact; instead, 
we can only show that the sequence J X \y\^dy is compact, from which we conclude a 
subsequence of S x that converges almost everywhere, thereby prove the existence. For 
the compactness proof in Section 4 we apply the compact Sobolev imbedding theorem, 
and don't need the Aubin-Lions Lemma or its generalized form of this lemma given by 
Roubicek [25], Simon [26], which plays a crucial role in the article [3]. 
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In the proof of regularity estimates, we differentiate equation (|1.2ip with respect to t. 
Thus a term like (\S x \ K )ttps appears, which cannot be absorbed by the a priori estimates 
with a weight \S X \ K . To overcome this difficulty, we derive a type of estimate (see (13. 6p ) 
with a reciprocal weight IS^I" 1 . This is possible due to the special structure of the model 
studied here. However the Allen-Calm model does not possesses such a structure, and 
our technique does not work for that model. 

The method of proof is limited to one space dimension, since for the a priori esti- 
mates it is crucial that the term | iS^c | J S f a . as in (jl.9p can be written in the divergence form 
^(\S X \S X ) X . In the higher dimensional case the corresponding term \V X S\A X S cannot be 
rewritten in this way, whence the multi-dimensional problem is still open. 



2 Existence of solutions to the modified problem 

In this section, we study the modified initial-boundary value problem and show that it has 
a Holder continuous classical solution. To formulate this problem, let \ £ Co°0^> [0, oo)) 
satisfy x(t)dt = 1. For k > 0, we set 



X n{t) := -X ( - 



and for S € L°°(Q Te ,M) we define 



( XK *S)(t,x)= X K(t-s)S(s,x)ds. (2.1) 
J 

The modified initial-boundary value problem consists of the equations 



-T lx = b, (2.2) 
T = D(e(u x )-e X n*S), (2.3) 
S t = cu\S x \ K S xx + c(T • e - ft(S))(\S x \ K - k), (2.4) 

which must hold in Qx e , and of the boundary and initial conditions 

u(t,x) = 0, (t, x) e (0, T e ) x dn, (2.5) 
S(t,x) = 0, (t,x) G (0,T e ) X dfl, (2.6) 
5(0, z ) = S (x), xefl. (2.7) 

To formulate an existence theorem for this problem we need some function spaces: For 
nonnegative integers m, n and a real number a £ (0, 1) we denote by C m+a (tt) the space 
of m— times differentiable functions on Q, whose m— th derivative is Holder continuous 
with exponent a. The space C a ' a ^ 2 (Q Te ) consists of all functions on Qt £ , which are 
Holder continuous in the parabolic distance 



d((t,x),(s,y)) := yj\t - s\ + \x - y\ 2 . 

C m ' n (Q T J and C m+a ' n+a / 2 (Q T J, respectively, are the spaces of functions, whose x— 
derivatives up to order m and i-derivatives up to order n belong to C(Q T J or to 
C a ' a ^ 2 (Q T J, respectively. 
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Theorem 2.1 Let u, k > 0, T e > 0. Suppose that the function b G C(Q Te ) has the 
derivative bt G C(Q Te ) and that the initial data Sq G C 2+a (Q) satisfy So\gn = So jX \an = 
So,xx\dfi = 0. Then there is a solution 

(u,T,S) G C 2 'HQ Te ) x C 1 ' 1 ^) x C 2+a > 1+a ' 2 (Q T J 

of the modified initial-boundary value problem \2. 2]) - (2.1). This solution satisfies St x £ 



L 2 (Q Te ) and 

max 1 5 1 < max | So | . (2-8) 

Proof. In [3] it is shown that the unique solution to the linear elliptic problem (|2.2j) - 
(|2.3p . with (|2.5p and given 5, is given by 

u(t,x) = u*n(xK*S)(t,y)dy-^^^( X K*S)(t,y)dy\+w(t,x), (2.9) 

De* f d 

T(t,x) = D(e* -e)( XK *S)(t,x)-- ( Xk * S)(t,y)dy + a(t,x), (2.10) 

d- a J a 

where u* G R 3 ,e* G <S 3 are suitable constants only depending on e and D, and where for 
every t G [0,T e ] the function (w(t),o~(t)) : f2 — >■ M 3 x <S 3 is the solution to the boundary 
value problem 

-o- lx {t) = b(t), 

a(t) = De(w x {t)), 
w(t)\dn = 0. 

Since by assumption b and bt belong to C(Qt £ ), it follows that (w,a) G C 2,1 (QT e ) X 
C^'HOtJ- We insert (|2.10p into ([2~2]) and obtain the equation 

S t = ai(S x )S xx + a 2 (t,x,S,S x ,XK.* S i~j~ J (Xk * S)(t,y)dy\ ( 2 -H) 
in Qt b , where 

«i(p) = 

and 

a2(t, SjP, r, s) = c (e ■ D(e* — e)r — e ■ De*s + e ■ a(t, x) — 4>'(S)j (\p\ K — n). 

Then ()2.11|) . (|2.6p and (|2.7p form an initial-boundary value problem with nonlocal terms, 
which is equivalent to the problem (|2.2|) - (|2.7|) . We can apply [221 Theorem 2.9, p. 23], 
with a slight modification, to (|2.1ip and conclude the existence of classical solution to 
the modified problem and estimate (I2.8P holds. For the details, we refer to [3]. 
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3 A priori estimates 



In this section we establish a-priori estimates for solutions of the modified problem, which 
are uniform with respect to k € (0, 1]. We remark that the estimates in Lemma 13. II and 
some in Corollary C3.2, though stated in the one-dimensional case, can be generalized 
to higher space dimensions. 

In what follows we assume that 

0<k<1, (3.1) 

since we consider the limit k — > 0. The L 2 (f2)-norm is denoted by || • ||, and the letter C 
stands for varies positive constants independent of k, while may depend on v. 

We start by constructing a family of approximate solutions to the modified problem. 
To this end let T e be a fixed positive number and choose for every k a function Sq € 
C£°(0) such that 

ll^o _ ^o|lHi(n)ni? 2 (Q) °> K -> °> ( 3 - 2 ) 

where So S Hq(£1) n i^ 2 (17) are the initial data given in Theorem ll.2i We insert for So in 
(12. 7p the function Sq and choose for b in (I2.2p the function given in Theorem 1.1. These 
functions satisfy the assumptions of Theorem 12. 1\ hence there is a solution (u K ,T K , S K ) 
of the modified problem (12.2p - (12.7j) . which exists in Qr e - The inequality (12.81) and 
Sobolev's imbedding theorem yield for this solution 

sup \\S K \\ L °o( Q ) < sup \\l°°(0) < C- (3-3) 

0<K<1 0<K<1 

Remembering that a in dZHDJI belongs to C^^QtJ, we conclude from (13. 3|) that also 

max|c(T K -e-^{S K ))\ < C. (3.4) 



Lemma 3.1 There hold for any t E [0,T e ] 

\S«(t)f + cv f [ \S«\ K \SZ x \ 2 dxdr < C, (3.5) 
Jo Jn 

ft r |5«i 2 



, -dxdr < C. (3.6) 
o Jn \ D x Ik 

It is worth a remark on the estimate (|3,6p . 
Remark 1. The reciprocal weight estimate (|3.6|) .' It is obtained by multiplying equation 
(|2.4p by the reciprocal of the weight function \S£\ K in the leading term. This is due to the 
special structure of our model which makes it possible to get the regularity results in this 
article. In contract, the model of the Allen-Cahn type with the mobility depending on the 
gradient of unknown does not possesses this feature, for instance, St — \S x \ K S xx +f(S) = 0. 
Suppose we obtain the L°°(QT e )-norm of S. If we multiply it by S^IS^I" 1 and integrate 
it with respect to x, we then see that the term f(S)St\S x \~ 1 dx cannot be absorbed by 
the left hand side. Consequently, it is impossible to obtain an estimate of the form (]3.6p 
for the Allen-Cahn type equation. 
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Proof. Invoking S^ x € L 2 (QtJ, by Theorem 12. l\ which yields that for almost all t 



Making use of this relation and estimate (|3.4p . multiplying (|2.4p by — and integrating 
it with respect to x, and taking the boundary condition (12.6f) into account, we obtain 
that for almost all t 

\j t \\SZ\\ 2 +cvj^ \S%\ K \SZ x \ 2 dx = J c(i,'(S K ) - T K ■ e)(\S£\ K - n)S* xx dx 

< C (\S X \ K + n)\S xx \dx = C / liS'aj \k \S x |k IS^cte + C / ^S^ldx 

7n Jn 

. CV f | y-yrei I nK 1 2 7 II II 2 ^ /* / 1 ok I \2 J /~1 

< — / lOj. UPzz I uX + — KlPzzll H / (.PasUJ «X + G^. 
4 Jf! 4 ^ ./fi 



(3.7) 



Splitting the second term on the left hand side of (|3.7p into two equal terms and sub- 
tracting the term ^ f Q \S x \ K \S xx \ 2 dx and ^k||5 , ^ x || 2 on both sides of this inequality, and 
using Gronwall's Lemma we derive (|3.5[) from the resulting estimate, noting also (|3.2p 
and k < \S X \ K . 

To derive (|3.6p . we multiply (|2.4p by 5f jS^I" 1 and integrate the resulting equation 
with respect to x to get 



= / Wr dx ~ c I (vSZ x -i>s)S?dx + c [ {^-Sfdx 
Jn Px k Jn Jn \°x l« 

(Sf) 2 



Qk\ 
°x \k 



dx + h+ J 2 . 



(3.8) 



Invoking the formula tps = — ^ • e + ip'(S) and the boundary conditions, and using 
integration by parts we have 



cj [ (^\S*\ 2 + i>(S K ))dx-c f T«-eS?dx 



(3.9) 



To deal with the term J, we use ([3.4p to get 



\J\ 



1 

T ' £ [ S x | ^ y C&E 

I Orel 2 



K\2 



£ C (115^11 + 1)2 



(gf) 

n I'-'x Ire 



da; 



I <?«l 

'-'a! « 



< 



1 /■ (5, 



re\2 



H Pas 



Orel 
°i re 



dx + C. 



(3.10) 
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Here we used the estimate (|3.5p and the Cauchy-Schwarz and Young inequalities. 
For 12, we make use of equation (|2.4p and write 



:J ^Ucu\S%S^ - <*I> S {\S% - K))dx 



in \ \S£U 
By definition, one has \S£\ K > K which implies rm - < 1. So 



dx. 



(3.11) 



\h\ < C f (K\S^\ + \S^\ K + K)dx 

Jn 

< cHI^I| 2 + a(||^|| + l). (3.12) 
With the help of f|3.5jl . (|3.9p - (|3.12p . we integrate (|3.8p with respect to t, then obtain 

i i* r (sn 2 




o Jn 



dx + c/ (^|5^| 2 + V(^ K )) ^ < C 



(3.13) 



which implies (|3.6p . Thus we complete the proof of this lemma. 
Furthermore, we obtain 

Corollary 3.2 There holds for any t € [0, T e ] 

' (\SZ\ K \S£ x \)UxdT < c, 




io Jn 

ft 




o Jn 



(\S«S2 x \)zdxdT < C, 



\y\ndy 



w 1[ 3(n) 



dr < C, 



\y\ndy 



dr < C, 



L°°{n) 



Ill x\ x H^f (o^x^n)) 



t 8 
||5«||3 



L°°(n) 



< c, 

dr < C. 



(3.14) 
(3.15) 

(3.16) 

(3.17) 
(3.18) 

(3.19) 



Proof. For some 2 > p > 1 we choose g, g' such that 

2 11, 
q = ~, - + -7 = 1. 
P 9 9 



10 



By Holder's inequality, we have 

f [ {\s«\ K \s« x \y dxdr 

Jo Jtt 



ft 

'o Jn 




< 



< 



(\SZ\ K )*((\S«\ K )*\S«J)dxdT 

ft (|^U)^ dxdrY ( f ! (\S2\ K )V \S^rdxdr) 
Jo Jn J \Jo Jn J 



2-p 

f J (\SZ\ K )& dxdA 2 (J f^\S K x \ K \S^\ 2 dxdry . (3.20) 



Estimate (|3.5p implies that if p satisfies < 2, i.e. p < |, then the right hand side of 
(|3.20p is bounded. This yields estimate (|3.14p . By definition of \y\ K , 

,■2 



\vU - \v\ 



\yU + \y\ 



Since \y\ K + \y\ > k, we have 



Hence 



K 2 , K 2 



< — = k. (3.21) 



\y\ K + \y\ k 
o < \y\ R - \y\ < k 



Letting y = S£ yields 

\^x^xx\ = l^xW^xxt — \\^x\k ~ \Sx \)\Sxx\ — K \Sxx\i 

and (|3.15|) follows from (|3.14|) and estimate (|3.5|) . 
Next we are going to prove (|3.16p . Writing 

\S2\kSZ, = ( [ X \y\ K dy] , (3.22) 



o 



and invoking that the primitive of \y\ K is equal to 

- (yy/y 2 + k 2 + K 2 log (y + \/y 2 + K 2 

which, thanks to logx < x — 1 for all x > 0, is bounded by C(y 2 + 1), we then show 
easily that 

I ( \y\ K dydx <C I ' (\S*\ 2 + l)dx <C. 
Jn Jo Jn 

To apply the Poincare inequality of the form 

11/ - /IUp(fi) < C\\fx\\Lp(n) 
where / := -py f Q f(x)dx, we choose 

4 fSS 
P = g, / = J \y\*dy, 
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and obtain 



\y\ndy 



L3(n) 



dr 







4 

3 ft 




<- c f 




dr + C 


/ \y\ K dy 


Jo 




L3(fi) J0 


Jo 



d.T 



L3(fi) 



< C I m\ K S^\\» L ^dr + C I ldr, 



which implies, by (|3.14|) . that 



f 


r s x 


\y\ndy 


10 


Jo 





dr <: C. 



(3.23) 



(3.24) 



L3(fi) 

Hence (I3TT6D follows, and we get / Q S " \y\ K dy € Li (0, T e ; VF 1 '! (0)). Making use of the 
Sobolev embedding theorem, we get (|3.17p . 

It remains to prove estimate (|3.19|) . since (|3.18|) is equivalent to (|3.19|) . We rewrite 
Jo* \y\ndy as 



\ K dy 



\y\dy+ / (\y\ K - \y\)dy 



\y\y 



\s%\s% + 



Thus 



\y\dy 



o bU + bl 
o bU + bl 



|yU<% 



K 2 S K 



(3.25) 



(3.26) 



By (|3.2ip and the Young inequality we obtain from (|3.5p and the assumption that < 1 
that 



\Sx\k + Pal 



< 1^1, thus 

3 

< I J (k 2 + k\S* x \ 2 ) dxdrj < C. 



(3.27) 



Combination with (|3.16p and (|3.26j) yields 



\m\s« x)x \\ LkQTe) <c 



\y\ndy 



o 



+ C\\kS^ x \\ a <C. (3.28) 



It is clear that \S^\S^ < C J n \S£\ dx < C. Applying again the Poincare inequality 
to the function / = \S£\S£, we arrive at 



II \s K x \s K x \ 



4 <C. 

L?(Q Te ) - 
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Hence this, combined with f)3.28p . implies that 

1 1 I QK I QKj II ^ /~Y 

II I XI a: "L^^TejW 1 '^^)) ~~ 
one concludes by using the Sobolev embedding theorem that 

1 1 | QK | QK || ^ S~1 

IN x u L?(0,T e ;L°°(n)) - ' 

which is 

\\Sx\\ « , „ < c. 

This completes the proof of this corollary. 
Lemma 3.3 There hold for any t 6 [0, T e ] 

||^ K (t)|| 2 + c^ J J \S^\ K \S^\ 2 dxdT < C, (3.29) 
\\\S2\A(t)\\ < C. (3.30) 

Proof. Suppose that estimate (|3.29|) is true, from equation (|2.4|) and the estimate 
we can easily get (|3.30p . So it is enough to prove (|3.29p . 
Differentiating (|2.4p formally with respect to t yields 



S£ = cv(\SZ\ K SZ t ) x + c((T K -e-ft(S K ))(\SZ\ K - K )) t . (3.31) 
Multiplying (|3.3ip by S* and integrating it, by integration by parts, we obtain 

= \jJ S t\\ 2 + ^f \S2\ K \S£ t \ 2 dx + c f (ip«-E-$(S K ))(\S2\ K -Kj) S?dx 

= ~ira 2 + a/y |5«U|5^| 2 dar + Ji. (3.32) 

It is not difficult to carry out a rigorous justification of (|3.32p with the help of difference 
quotient, we omit the details. Computation gives 



Ji = C y^((T K - £ --^(^)) t (|5«U- K ) + (T K . e --^(5 K ))(|^U) i ) < S i K dx 

= J11 + J12. (3.33) 
By the formula of T, we have 

I Jul < C f (\S?\ 2 + \S?\)(\SZ\ K + K)dx 
Jn 

< C(||S£|| ioo(n) + l) / (\S?\ 2 + l)dx 

Jn 

< C(\\S£\\ Loa{Q) + l)(\\S?f + l). (3.34) 
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To handle J12, we make use of estimate (|3.6p and \y\ < \y\ K . 



12 



< C 



< c 



I QK CK QK I 
\ D X D Xt°t I 

Pa; U 



dx = C 



I Q K I I C K 
Pa;°a:tl Pt 



\ J x\K \ j x\k 



-dx 



z Pa; l« 

Thus it follows from (13732]) - (13735]) that 



(3.35) 



cv 



< T \S*\ K \S^dx + C u ^-dx + C(||^|| Ji< » (n) + l)(||5fr + l).(3.36) 

From equation (|2.4p and assumption So € H 2 (£l) we compute the initial data 

||S?|t=o|| < C{\\\S 0x \ k S Qxx \\ + \\\S Qx \ k + k\\) 

< c(\\\s 0x \ II^Oxxl! + II^Oill + 1) 

< C((\\S 0x \\ H i (n) + l)\\Soxx\\ + \\S 0x \\ + l)<C. (3.37) 

Thus St \t=o £ L 2 (fi). Next we use the Gronwall inequality of the form: 

Lemma 3.4 For measurable functions y,A,B defined on [0,T e ], such that y > and 
A : BeL x (Q,T e ),if 

y'(t)<A(t)y(t) + B(t), 
then ^ ^ 

y(t) < y(0) exp (J Air )dr "\ + J B(s) exp (J A(t )dr "\ ds. 

Defining 

y(t) = \\S?(t)\\ 2 , A(t) = C(\\S£\\ LOO{CI) + 1), B(t) = (7(11^11^(0) + 1) + C V 



K\2 



\S£\ 



-dx, 



where A, B are integrable over [0,T e ] by Lemma 13. II and Corollary 13.21 we derive from 
(13736]) and (137371) that 

WSmf + w f I \S*\ K \SZ t \ 2 dxdT<C\\S?(0)\\ 2 + C<C. (3.38) 
Jo Jn 

Thus the proof of this lemma is complete. 



Corollary 3.5 The function \y\ K dy belongs to H 1 (Qt e ), and the estimates hold 



\y\Zdy 



\%dy 



<C. 



L 2 (Q Te ) 



<c. 



L 2 {0,T e ;H^(Q)) 



(3.39) 
(3.40) 
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Proof. Calculating yields 



X 1 

y\ldy 



recalling (|3.29p . we obtain (|3.39|) . Similarly, 



combining with (|3.5p gives 



x 1 

y\*dy 



OK 

°X 1 



I QK I 2 QK 



I CK I 2 CK 







< c. 



(3.41) 



(3.42) 



(3.43) 



L 2 (Q Te ) 



have 



Finally, Noting \ Jq x \y\Zdy\ < Cmax |m, 2 1 for some large constant M > 0, 

/ \S£\ 3 dx 
Jn 



we 



LHQ Te ) 



< c + c 



< C + C\\S£\\Loo(fy I \Sx\ 2 dx. 



Thus by (|3.19j) in Corollary 13.21 there holds 



^x 1 

\y\Zdy 







oIt < C + C \\S: 



're llL°°(n) 



dr < C. 



(3.44) 



(3.45) 



L 2 (Q Te ) 

Then ()3.40p follows from (|3.43p and ()3.45p . The proof of the lemma is complete. 



4 Existence/regularity of solutions to the phase field model 

We shall make use of the a priori estimates established in the previous section to study 
the convergence of {u K ,T K , S K ) as k — >• 0. In this section we will show that there is a 
subsequence, which converges to a weak solution of the initial-boundary value problem 
(]1.7p - (]1.12p . thus we prove the existence of weak solutions; then we shall investigate 
the regularity of solutions. 

Existence. It follows from Lemmas 13.11 and 13.31 that 

\\S k \\ hHQtJ <C, (4.1) 

for a constant C independent of k. Hence, we can select a sequence n n — > and a function 

5 € H 1 {Qt<,), such that the sequence S Kn , which we again denote by S K , satisfies 

\\S K - S || L 2 (Q Te ) ~^ °> S^^S X , Sf^St, (4.2) 

where the weak convergence is in L 2 (Qt e ) ■ 

Since the equation (II. 9p is nonlinear, the weak convergence of S£ is not enough to 
prove that the limit function solves this equation. In the following lemma we therefore 
show that converges pointwise almost everywhere: 
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Lemma 4.1 There exists a subsequence of S£, we still denote it by S£, such that 
f S: - 

J \y\%dy ->• / \y\ 2 dy a.e. in Q Te , (4.3) 

a.e. in Q Te , (4.4) 





l 

\y\ 2 dy 




Jo 




QK 
°X 




QK 1 


\S x \i 


rS x 


1 

y\ 2 dy, 


Jo 



weakly in L 2 {Qx e ), (4-5) 

/ \y\ 2 dy — > / \y\zdy, strongly in L 2 (Qt £ ), (4-6) 
Jo Jo 

as k — y 0. 

The proof is based on the following result: 

Lemma 4.2 Let (0,T e ) x £1 be an open set in M + x W 1 . Suppose functions g n ,g o,re in 
L 9 ((0,T e ) x D,) for any given 1 < q < oo, which satisfy 

\\9n\\Li((o,T e )xn) <C, g„^ g almost everywhere in (0,T e ) x tt. 

Then g n converges to g weakly in L 9 ((0,T e ) x fi). 

Since we have stronger a priori estimates than those in [3] where it is assumed that the 
initial data is in H 1 ^), we don't need any more a compactness lemma of Aubin-Lions 
type or its generalized version (see e.g. Simon |26] and Roubfcek |25j ) which in [3 1 plays 
a crucial role in the proof of the existence of weak solutions with H 1 ^) initial data. A 
proof of Lemma 14.21 can be found e.g. in the book by Lions |23} p. 12]. 

Proof of Lemma \4-l\ Since the estimates in Lemma 13.11 and Corollary 13.51 imply that 
the sequence L x \y\ndy is uniformly bounded in H (Qt b ) f° r K € (0, 1]. By the Sobolev 

OK I 

imbedding theorem, we assert that f Q x \y\Zdy is compact in L 2 {Qt e ) = L 2 (0,T e ; L 2 (Cl)). 
Thus there is a subsequence, still denoted by J x \y\K.dy, which converges strongly in 

L 2 (Qt £ ) to a limit function G € L 2 (Qt £ )- Next we prove that the sequence Jq x \y\^dy 
converges to G in L 2 (QtJ- Write 

°r J f°X 1 f°X I 1 







\y\ 2 dy = J {ylndy + J [\y\2 - \y\% )dy = L 1 +J 2 . 
It is easy to compute that 

n ^ i ,5 I 1 1 \y\ K - \y\ 
< \y\n -\y\ 2 = 



,i ,i 
\y\K + \y\ 2 

K 2 

(\y\l + |y|')(|yU + |y| 

- 2 



< -r TT = ^- (4.7) 



K2 



Thus I2 can be estimated as 

\\h\\tf(Q Te ) < \\K 2 S x \\ L 2 { Q Te) < CK2\\S£\\ LaO ( 0)Te . L a(ciy < Ck2 -> 0. 
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Therefore, f^ x \y\?dy — > lim re _>o h = G strongly in L 2 {QtJ). 

Consequently, from this sequence J x \y\2dy we can select another subsequence, de- 
noted in the same way, which converges almost everywhere in Qt e - Using that the 
mapping y 1— > f(y) := \y\2dy has a continuous inverse / _1 : R — > R, we infer that also 
the sequence S x = f^ 1 ( jF x \y\^dy ) converges pointwise almost everywhere to f~ l (G) 



Jo 

in Qt s - From the uniqueness of the weak limit we conclude that f^ 1 (G) = S x almost 
everywhere in Qt £ - Thus we prove (|4.6p . 

To prove (I4.5P we note that the estimate \S£\ K < \S£\ + k and the inequality (I4.ip 
together imply that the sequence \S£\ K is uniformly bounded in L 2 (Qt e )- Thus, (j4.5j) is 
a consequence of (14. 4h and Lemma 14.21 

Proof of Theorem \1.2c Define the functions u,T by inserting S into (|2.9[) and (|2.10|) . 
respectively, where S is the limit function of the sequence S K . We shall prove that 
(u, T, S) is a weak solution of problem (|1.7|) - (|1.12p . 

Recalling (12. 8p we have S* £ L°°(Qj i e ). From this relation, from the definition of u, T 
we immediately see that «, T satisfy (|1.14p and (|1.15|) . respectively. Observe next that 
ll^llz^o.Tej.H^fi)) — G, by Lemma I3TT1 and Sobolev's embedding theorem. This implies 
S € L°°(0, T e ; Hq(Q) ), since we can select a subsequence of S K which converges weakly 
to S in this space. Thus, S satisfies (|1.16p . 

Noting that from (pTT]) and (14^21 



||x K *5 K -5|| L2(QTe) < \\XK*(S K -S)\\ LHQTe) + \\(S-x K *S)\\ L 2 iQTe) 

< \\(S- X .*S)\\ L2{QTe) + \\S«-S\\ L2(QTe) ^0, (4.8) 

for k — > 0, we conclude easily that the function (u, T) defined in this way satisfy weakly 
equation (|1.7p - (jl.gp . By definition, if the relation (|1.17p holds, then the proof of the 
existence of weak solutions is complete. To verify (|1.17p we use that by construction S K 
solves (12. 3j) . Now we multiply equation (|2.3j) by a test function ip G Cq°((— 00, T e ) x O) 
and integrate the resulting equation over Qt e , then obtain 

= (S?,v) QTe +(-cv\S K J K S K xx + F K (\SZ\ K -K),<p) QT 
= -(S$,<p(0))n - (S K ,vt)Q Te + (cu f ' \y\ K dy,<p x 



Q 



+ (T k (\S*\ k -k)^) Qt , 

where T n = —c{T K ■ e — vp'(S K )). Equation (|1.17p follows from this relation if we show 
that 

(SS,if(0)h (SoMtyfo, (4.9) 
(S^tW (S,ipt) QT , (4.10) 



y\ K dy, ip x ) — I -\S x \S x ,<p x 1 . (4.11) 



/Q Te V 2 /Qt, 



(^|^U,(^) QTe -> (^|5 x .U) QTe , (4.12) 
( K 7- K ,V;) QTe -> 0, (4.13) 
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for k — y 0. Now, the relation (14. 9p follows from (j3.2j) . and the relation (|4.10|) is a conse- 
quence of (|4.2j) . By 114.4 jl and (13.301) from which it is easy to get || \y\ K dy\\i2(Q T j < C, 
using again Lemma 14.11 one has (|4.1ip . Convergence f|4. 13|) follows from (|3.4p easily. 
To verify KW\ we note that gSJ), |U} , flUJ) , and the definition of J" K yield 

||^|^U|| L2(Qre) < C, (4.14) 

J rK |5'xU ~~ ^ -^"l^l) almost everywhere. (4-15) 

Then by Lemma |4.H 

-TK I I ^ 77| O I 

weakly in L 2 (Qr e ), which implies (|4.12p . Consequently (|1.17p holds. 

Regularity. Since So G H 2 (£2), we can obtain more regular solutions. By the estimate 
l|5't c |lL oo (o,T e ;L 2 (r2)) — C> we see that the relation St G L°°(0, T e ; L 2 (S7)) is true. Then by 
the theory of elliptic systems, we obtain (|1.19p . 

To prove (ll.20p . we recall the definition of weak solutions. From (|1.17p it follows that 



\S X \S X , <Px)qt. I — c 



(^(T-e-^(S))\s x \,<p) ( 



+ \(S,<p t ) QTi .+(S M0))a\ 



< C||5'a;||L° o (0,T e ;L 2 (n))ll9 : 'llL 1 (0,T e ;L 2 (n)) + I (<% j ^)Qt s I 

< C|MU*(0,T e ;L 2 (fi)) + ll»S't||L°°(0,T e ;L 2 (n))ll^llL 1 (0,T e ;L 2 (n)) 

< C\\(p\\ L l(p )Te . L 2(p^, (4-16) 

here, we used the estimates ||<S'x||i°o(o,T e ;L 2 (Q)) < C an d ||<St ||L°°(o,T e ;L 2 (f2)) < C- The right 
hand side of (|4"T6l) holds for all 99 G L*(0, T e ; L 2 (ft)), whence 

sup HGSxls^WII = sup \((\s x \s x ) x ,(f) QTe \ 

°^ T e IMI.tl(0,T e ;£2(C2))<l 

= sup ^ 1(1^15^, ^)qt 6 I 

< C. (4.17) 

Thus, (\S X \S X ) X G L°°(0,T e ;L 2 (J))). 

Furthermore, from the Poincare inequality and the estimate ||5' :r ||ioo(o i T e ;L 2 (f2)) — C 
we obtain (S^S^ G L°°(0, T e ; H 1 ^)), from which one asserts by the Sobolev imbedding 
theorem that \S X \S X G L°°(Q Te ), hence S x G L°°(Q Te ). And the proof of Theorem [L2l is 
complete. 
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